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Abstract 

Introduction: Phosphatidylinositol-4,5-bisphosphate {PIP2) is a cofactor necessary for the activity of KCNQl channels. Some 
Long QT mutations of KCNQl, including R243H, R539W and R555C have been shown to decrease KCNQl interaction with 
PIP2. A previous study suggested that R539W is paradoxically less sensitive to intracellular magnesium inhibition than the 
WT channel, despite a decreased interaction with PIP2. In the present study, we confirm this peculiar behavior of R539W and 
suggest a molecular mechanism underlying it. 

Methods and Results: COS-7 cells were transfected with WT or mutated KCNEl -KCNQl channel, and patch-clamp 
recordings were performed in giant-patch, permeabilized-patch or ruptured-patch configuration. Similar to other channels 
with a decreased PIP2 affinity, we observed that the R243H and R555C mutations lead to an accelerated current rundown 
when membrane PIP2 levels are decreasing. As opposed to R243H and R555C mutants, R539W is not more but rather less 
sensitive to PIP2 decrease than the WT channel. A molecular model of a fragment of the KCNQl C-terminus and the 
membrane bilayer suggested that a potential novel interaction of R539W with cholesterol stabilizes the channel opening 
and hence prevents rundown upon PIP2 depletion. We then carried out the same rundown experiments under cholesterol 
depletion and observed an accelerated R539W rundown that is consistent with this model. 

Conclusions: We show for the first time that a mutation may shift the channel interaction with PIP2 to a preference for 
cholesterol. This de novo interaction wanes the sensitivity to PIP2 variations, showing that a mutated channel with a 
decreased affinity to PIP2 could paradoxically present a slowed current rundown compared to the WT channel. This suggests 
that caution is required when using measurements of current rundown as an indicator to compare WT and mutant channel 
PIP2 sensitivity. 
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Introduction 

Cholesterol regulates several ion channels, and changes in 
membrane cholesterol levels provoke various effects depending on 
the channel type (reviewed in [1,2]). Regarding inwardly rectifying 
channels (Kir), cholesterol effects have been studied in detail. 
Effects of cholesterol increase or depletion vary among families 
and even among channels of the same family. Similarly to Kir 
channels, functional effects of cholesterol on voltage-gated (Kv) 



channel activity are highly variable: an increase in Kv2.1 current 
in DrosophUa neurons [3] and a suppression of the same current in 
mammalian pancreatic |3-cells [4]. Adding more complexity, 
effects on the cardiac potassium channel KCNQ^l are variable 
depending on the drug used to decrease cholesterol levels [5,6], 
More specifically, Probucol - known as a cholesterol depleting 
agent - is able to decrease the coexpressed KCNEl /KCNQ^l 
current amplitude [5] without decreasing cholesterol levels in 
CHO-Kl [6]. Simvastatin and triparanol are more specific since 
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their main efTects (activation kinetics) are similar and correlated to 
a cholesterol decrease. Clearly, cholesterol effects on KCNQ_I 
channels remain to be studied in more detail. 

Phosphatidylinositol-4,5-bisphosphate (PIP^) is abundant in 
cholesterol-rich membrane domains [7]. The mechanisms by 
which PIP2 regulates several channels, including KCNQl, have 
been studied in much greater detail than the mechanism of their 
regulation by cholesterol [8-16]. PIPj is a minor acidic membrane 
lipid found primarily in the inner leaflet of the plasma membrane. 
It has been shown to be a necessary cofactor for a wide variety of 
ion channels and transporters [17]. In Kir and Kv channels, 
several stimuli impact channel activity by decreasing available 
PIP2 or modulating channel-PIP2 interactions [8,9,18-20]. Con- 
sistent with PIPj regulating channel activities, mutations that 
impair channel-PIPj interactions play a crucial role in channelo- 
pathies. Regarding KCNQ^l , we have shown that type 1 long QT 
(LQTl) syndrome can be associated with a decrease in KCNQ^l- 
PIP2 interactions provoked by mutations in the S4— S5 linker 
(R243H) and in the C-terminal domain CTD (R539W and 
R555C) [10]. This has since been confirmed for mutations 
R539W and R555C [11,12]. 

In previous studies, PIP2 sensitivity of the KCNQ_1 mutants 
R243H, R539W and R555C was determined by using the soluble 
short acyl chains analog diC8-PIP2 [10]. Since then, PlPg- 
regulated Kvll.l channels have been shown to be insensitive to 
diC8-PIP2 [2 1] , suggesting that channel affmity for PIP2 and diC8- 
PIP2 might be quite different. Moreover, the R539W behavior 
seems paradoxical since, in COS-7 cells, this mutant was shown to 
be less sensitive than WT to intracellular Mg^"^, which decreases 
PIP2 availability through masking its negative charges [10]. In 
general, any mutation that reduces channel affinity to PIP2 
provokes an increased inhibition by Mg^^ [18]. Hence, we set out 
to use other approaches to further study the PIP2 affinity of 
KCNC)1 mutant channels in the presence of KCNEl. We 
measured the kinetics of Mg^"""-, wortmannin- and Ci-VSP- 
induced current rundown and confirmed a decreased PIP2 aflOnity 
for two of the three mutants. To our surprise, for the third mutant 
(R539W) the current is running down more slowly when a\'ailable 
PIP2 is decreased. We therefore hypothesized that the KCNQ^l- 
R539W mutation leads to the stabilization of an open-pore 
conformation shortcutting the PIP2 effect on the concerted 
opening. Using both modeling and functional analyses, we show 
for the first time that a mutation in the CTD domain shifts the 
channel interaction with PIP2 to a preference for cholesterol. 

In conclusion, our study shows that cholesterol affects the 
rundown of R539W but not of WT KCNE1-KCNC)1 channels, 
suggesting that cholesterol specifically stabilizes opening of 
R539W channels and decreases their need for PIP2 to be open. 
Estimating a channel affinity to PIP2 from the kinetics of the 
current rundown when PIP2 decreases is an approach that has 
been commonly used. This study indicates that such an approach 
is not always relevant, since the mutation that disrupts the 
interaction with PIP2 may stabilize the channel open state through 
another mechanism. 

Methods 

Cell culture and transfection 

The African green monkey kidney-derived cell line COS-7 was 
obtained from the American Type Culture Collection (RockviUe, 
MD, USA) and cultured in Dulbecco's modified Eagle's medium 
(GIBCO, Paisley, Scotland). Cells were transfected with plasmids 
(2 |ig per 35 mm dish) complexed with Fugene-6 (Roche 
Molecular Biochemical) for giant-patch and ruptured-patch 



experiments, or Jet-PEI (Polyplus-Transfection) for permeabi- 
Uzed-patch experiments, according to the standard protocol 
recommended by the manufacturers. We used a construct of 
hKCNQ,! fused to hKCNEl (KCNEl -KCNQ,1) to prevent 
variability caused by a variable KCNEl /KCNQ,1 expression 
ratio [10]. This type of construct is justified by the fact that the 
cardiac channel is generated by the assembly of 4 KCNQ_1 and up 
to 4 KCNEl subunits [22]. For giant-patch experiments, the 
relative DNA composition was 80% pCDNAS.l plasmid contain- 
ing the human WT or mutated KCNEl -KCNQ^l concatemer 
[10,23] and 20% pEGFP coding for the green fluorescent protein 
(Clontech). For permeabilized-patch experiments testing the effects 
of osmolarity, relative DNA composition was 40% pCDNA3.1 
KCNEl -KCNQ,1 and 60% pEGFP. For permeabilized-patch 
experiments testing the PKA-dependent channel upregulation, 
relative DNA composition was 20% pCDNA3. 1 -KCNEl - 
KCNQ^l, 40% pEGFP and 40% pCDNA3-yotiao (a kind gift of 
Dr Robert S. Kass, Department of Pharmacology, College of 
Physicians & Surgeons, Columbia University, New York, NY, 
USA). For ruptured-patch experiments testing the effects of Ci- 
VSP, relative DNA composition was 20% pCDNA3. 1 -KCNEl - 
KCNQ^l, 70% pEGFP and 10% pRFP-Cl-Ci-VSP (a kind gift of 
Dr Dominik Oliver, Institute of Physiology and Pathophysiology, 
PhUipps-Universitat Marburg, Marburg, Germany). 

Electrophyslology 

Forty-eight to seventy-two hours post-transfection, COS-7 cells 
were mounted on the stage of an inverted microscope' and 
constantiy perfused at a rate of 2 mL/min. Experiments were 
performed at room temperature for giant-patch and ruptured- 
patch configurations, or 35.0±1.0°C for permeabilized-patch 
configuration. Stimulation, data recording, and analysis were 
performed either by Acrjuisl (Bio-logic Science Instruments, 
Claix, France) through an analog-to-digital converter (Tecmar 
TMIOO Labmastcr; Scientific Solution, Solon, OH, USA), or by 
pClamplO.4 through an analog-to-digital converter (Axon Digi- 
data 1440, Molecular Devices, Sunnyvale, CA, USA). Electrodes 
were electrically connected to a patch-clamp amplifier (RK-400, 
Bio-logic Scienc:e Instruments, Claix, France). 

For giant-patch experiments, the procedure has been described 
previously [24]. An 'excision' pipette, filled with the standard 
solution, was connected to a 10 ml syringe to apply suction for 
excision. Pipettes were pulled from borosilicate glass capillaries 
(glass type 8250; Gamer Glass) on a vertical puller (P30; Sutter 
Instruments Co., Novato, CA, USA) and fire polished using a 
microforge (MF-83; Narishige, Japan) to reach 9 to 12 |im tip 
diameters for patch pipettes and around 15 |J,m for excision 
pipettes. Cells were continuously perfused with the standard 
solution. KCNEl -KCNQl currents were analyzed using a 
protocol consisting of depolarizing voltage steps of 1 s from a 
holding potential of -80 mV to +80 mV and then to -40 mV for 
0.5 s, every 5 or 2 s. A microperfusion system allowed local 
application and rapid change of the different experimental 
solutions. 

For permeabilized-patch experiments, micropipettes (tip resis- 
tance: 2—3 MOhms) were pulled from soda lime glass (Kimble; 
Vineland, New Jersey, USA). KCNEl -KCNQ,1 currents were 
investigated with a protocol consisting of depolarizing voltage steps 
of 4 s from a holding potential of — 80 mV to +80 mV and then to 
—40 mV for 1 s, every 7 s. To obtain the WT activation curve, 
the membrane potential was stepped, from a holding potential 
(-80 mV) to six voltage steps (-20 to 80 mV, with a 20 mV 
increment) and then stepped back to —40 mV, where tail currents 
are visible. To obtain the R539W activation curve, the membrane 
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potential was stepped, from a holding potential (—80 mV) to six 
voltage steps (20 to 120 mV, with a 20 mV increment) and then 
stepped back to —40 mV, where tail currents are visible. 
Activation curves were fitted by a Boltzmann equation. KCNEl- 
KCNQ_1 deactivation kinetics were obtained by a monoexponen- 
tial fit. 

For ruptured-patch experiments (Ci-VSP), the same pipettes as 
for permeabUized-patch experiments were used. KCNE1-KCNQ,1 
currents were investigated with a protocol consisting of depolar- 
izing voltage steps of 2 s from a holding potential of —80 mV to 
+80 mV and then to —40 mV for 1 s, every 8 s. 

Patch-clamp data are presented as mean ± SEM. Statistical 
significance of the observed effects was assessed by one-way 
ANOVA or Student's t-tests. Off-line analysis was performed using 
Acquis 1, Clampfit and Microsoft Excel programs. Microsoft 
Solver was used to fit data by a least-square algorithm. 

Solutions and drugs 

For giant-patch experiments, the cells were perfused with a 
standard solution containing (in mmol/L): 145 KCl, 10 HEPES, 1 
EGTA (pH 7.3 with KOH). The following solution (in mmol/L): 
145 K-gluconate, 10 HEPES, 1 EGTA, (pH 7.3 with KOH) was 
used to perfuse the cell during K"*" current measurements and to fill 
the patch pipette tip (the CI -containing solution was in contact 
with the Ag^/AgCl filament). 

For permeabilized-patch experiments, the pipette (intracellular) 
solution had the following composition (in mmol/L): 145 KCl, 10 
HEPES, 1 EGTA and 0.85 amphotericin B (pH 7.2 with KOH). 
The Tyrode superperfusion solution contained (in mmol/L): 145 
NaCl, 4 KCl, 1 CaCl2, 1 MgClj, 5 HEPES, 5 glucose (pH 7.4 
with NaOH). The locally applied extracellular control solution 
(334 mosmol/L) contained (in mmol/L): 145 NaCl, 4 KCl, 1 
CaCla, 1 MgCls, 5 HEPES, 5 glucose, 20 mannitol (pH 7.4 with 
NaOH). Hypoosmotic challenge (234 mosmol/L) was induced by 
a decrease of NaCl from 145 mmol/L to 86 mmol/L and an 
increase of mannitol from 20 mmol/L to 39 mmol/L. In 
hyperosmotic solution (434 mosmol/L) NaCl was increased to 
205 mmol/L and mannitol was decreased to 0.9 mmol/L. We 
looked at the effect of osmolarity on R539W and WT channels in 
the same set of experiments. The effects on the WT channel have 
been previously published and ser\'e as a positi\'e control [14]. For 
ruptured-patch experiments, the pipette (intracellular) solution had 
the following composition (in mmol/L): 145 KCl, 10 HEPES, 1 
EGTA and 1 MgClj (pH 7.2 with KOH). CeUs were perfused 
with Tyrode solution. 

Available PIP2 was decreased by three ways: in the permeabi- 
lized-patch configuration, PIP2 was decreased by application of 
10 |J,mol/L wortmannin [16,25,26]. In giant-patch experiments, 
available PIPg was decreased by direct application of 1 . 1 mmol/L 
free Mg^"*" on the intracellular side of giant patches [14]. In 
ruptured-patch configuration, channels were coexpressed with the 
voltage-dependent phosphatase from Ciona intestinalis, Ci-VSP, 
allowing PIP2 dephosphorylation [27,28]. 

In channel PKA-dependent phosphorylation experiments, cells 
were exposed to a solution containing (in nmol/L) 400 cpt-cAMP, 
10 forskolin and 0.2 okadaic acid. For cholesterol depletion 
experiments, two different approaches were used: one-hour pre- 
treatment of 2 mmol/L 2-hydroxypropyl-P-cyclodextrin in Tyr- 
ode solution at 37°C or twenty-four hour pre-treatment of 
10 |a,mol/L triparanol [6] in 1 ml of the cell culture medium. 
Frc<; magnesium activities in gluconate-containing solutions were 
calculated using software designed by G. L. Smith (University of 
Glasgow, Scotland), using stability constants [29]. 



Molecular Modeling 

Peptides. 3D models of the WT and mutated sequence 
fragments, Q,QARKPYDV-R539-DVIEQYSQG and QQARK- 
PYDV-W539-DVIEQYSQG, were calculated using Peplook [30]. 
Briefly, populations of 3D models are built from sequences using a 
boltzmanian stochastic procedure in which calculated populations 
are 3D models of lower energy. Models are all-atoms and their 
energies take into account all interactions between non-bound 
atoms with a cut-off of 20 A. Energy terms are van der waals 
(Lennard-Jones), electrostatic (coulomb with an exponentially 
variable dielectric constant from 1 to 80 on a distance of 2 to 
30 A) and two hydrophobicity terms, one for intramolecular 
interactions and the second for solvent effect. The first term uses 
atomic transfer energy (Etr) and atom distances, the second one 
uses atom accessible atom surface ar(;a and atom Etr in the 
corresponding s()lv(;nt, as previously described [30]. The best 
model is named the Prime and the next 98 models of low energy 
are also sorted to form the peptide population. The populations 
are clustered around lead models on the basis of backbone 
RMSd<l A. 

Lipids. 3D models of lipids were calculated by the same 
procedure as PepLook but using a dataset of rotation angles 
covering the 360° rotation possibilities by steps of 7.5°. 

Impala. The Prime models of each peptide and Upid were 
individually tested across the water/membrane slab named 
IMPALA to define their best position [31]. IMPALA is a 
continuous restraint potential mimicking membrane hydropho- 
bicity, charge density and fluidity properties with respect to water. 
The mass centre of each model was set at 101 different levels i.e. 
every angstrom from —50 to -1-50 from the membrane centre, 
across the slab. At each position, the restraint potential of 10,000 
different orientations (rotations in the x/y plane of the membrane) 
was calculated and used to select the best position. In the figures, 
colored grids indicate different membrane interfaces: pink, the 
water/membrane interface; purple, the polar head/acyl chain 
interface of lipids; yellow, the membrane centre. 

Molecule hydrophobicity and charge. To illustrate mole- 
cule hydrophobicity and charge, we used MHP (Molecular 
Hydrophobicity Potential) and MEP (Molecular Electrostatic 
Potential). MHP and MEP are three-dimensional plots of the 
hydrophobicity and electrostatic isopotential surface of a molecule, 
respectively [32]. In the MHP plots, surfaces joining afl points of 
±0.1 kcal are drawn; green for the hydrophilic surface, brown for 
the hydrophobic surface. In the MEP plots, surfaces joining aU 
points of ± 10 kcal are drawn; blue for the negatively charged 
surface, red for the positively charged surface. 

Results 

R539W KCNE1-KCNQ1 is poorly sensitive to PIP2 decrease 

To obtain a quantification of channel PIP2 sensitivity, addition 
of soluble diC8-PIP2 is commonly used [8]. In a previous study, we 
used this approach to determine the PIP2 sensitivity of three 
LQTl -associated KCNQ,1 channels containing the R243H, 
R539W or R555C mutation [10]. However, this approach may 
be limited because, for some channels, the maximum effect of 
short-chain PIP2 is less than that of PIP2, suggesting a lower 
affinity [33] and/or a contribution of the acyl chain to the 
channel-PIP2 interaction. As an extreme example, the Kvll.l 
channel is sensitive to PIP2 but not to diC8-PIP2 [21]. These new 
results prompted us to use alternative methods to confirm the 
decreased PIP2 sensitivity of R243H, R539W and R555C mutant 
channels. We used three different approaches: Extracellular 
wortmannin effects on whole-ceU currents, intraceUular magne- 
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sium effects on excised-patch currents and Ci-VSP-coexpression 
effects on whole-ceU currents [10,18,27,28,34]. 

In the permeabilized-patch configuration, we monitored the 
activity of WT and mutant KCNEl-KCNQ^l channel currents 
during application of 1 0 |imol/L wortmannin that blocks the PI4- 
kinase required for PIP2 synthesis. Current density measured at 
the end of the depolarizing step (+80 mV) was monitored every 7 s 
and normalized to the current value measured before wortmannin 
application (time 0). Wortmannin application led to a gradual 
decrease of WT KCNEl-KCNQ^l channel activity, caUed 
rundown (Figure lA and IB). This rundown is also significant 
for the R243H and R555C, but not for the R539W channel 
(Table SI in File SI). In order to compare the rundown kinetics in 
the different conditions, the relative current amplitude after 63 s in 
10 |J,mol/L wortmannin was calculated. Rundown was accelerated 
for R243H and R555C mutants as compared to WT (Figures IB 
and IC), confirming a decreased PIP2 sensitivity. Relative current 
amplitude after 63 s in 10 (imol/L wortmannin was 0.88±0.03 
(n=7) for WT channel, and 0.77±0.04 (n = 5; P<0.05) and 
0.76±0.06 (n = 5; P<0.05) for R243H and R555C mutants, 
respectively. Most importantly, despite the lower affmity of short- 
chain PIP2 for R539W, the decrease of PIP2 levels caused by 
wortmannin application had no effect on this mutant channel 
activity (Figures IB and IC). 



We tried to decrease PIP2 intracellular levels to a higher extent, 
in an attempt to decrease R539W channel activity. To do that, we 
used intracellular Mg^"*", which is known to mask PIP2 negative 
charges [35]. We recorded WT and mutant channels activities 
during 1 . 1 -mmol/L free Mg^"*" application using the excised-patch 
configuration (Figure 2, Table S2 in File SI). As previously shown, 
Mg^'"' application led to a gradual decrease of WT KCNEl- 
KCNQ,1 channel activity (Figure 2A and 2B; T = 18±3 s, n = 12). 
R243H and R555C currents decreased faster than that of WT, as 
during wortmannin application (Figure 2B and 2C; R243H, 
T = 6.3±1.3s, P<0.01 and R555C, t = 4.7±0.5s; n= 10-12, 
P< 0.001 as compared to WT). Conversely, the R539W channel 
current ran down more slowly than the WT (t = 60±22 s; n = 8; 
P<0.05). 

To confirm these results, we depleted PIP2 by coexpressing the 
channels with the voltage-dependent phosphatase Ci-VSP, which 
is known to dephosphorylate PIP2 under membrane depolariza- 
tion [27]. Ci-VSP activation led to a gradual decrease of WT 
KCNEl-KCNQl channel activity (Figure 3 A and 3B). As a 
negative control, the WT KCNEl-KCNQl channel was ex- 
pressed without Ci-VSP and no current rundown was observed 
(Figure 3, Table S3 in File SI). Tail current density measured at 
— 40 mV was monitored and normalized to the current value 
measured at time 0. The relative current amplitude after 64 s of 
Ci-VSP activity was calculated. Rundown was accelerated for 
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Figure 1. R539W is insensitive to wortmannin. A, representative permeabilized-patch current recordings of WT or mutant channels measured 
before or after a 63-s wortmannin application (10 |imol/L), and during the voltage-clamp protocol shown. Start-to-start interval = 7 s. B, relative 
current amplitude of WT or mutant channels measured at the end of the depolarizing step (4-80 mV), plotted against time. Current values are 
normalized to the current level measured before wortmannin application (time 0). These experiments were performed at 35°C in the permeabilized- 
patch configuration. In this configuration, it has been shown that there is no spontaneous current rundown [24]. C, mean relative current amplitude 
of WT or mutant channels measured after a 63-s wortmannin application (n = 5-7). *p<0.05, versus WT. 
doi:10.1371/journal.pone.0093255.g001 
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Figure 2. R539W is poorly sensitive to intracellular magnesium. A, representative giant-patch current recordings of WT or mutant channels 
measured after 5 and 25-s iVlg^* application (1.1 mmol/L free Mg^*), and during the voltage-clamp protocol shown. Start-to-start interval = 5 s. B, 
relative tail-current amplitude (measured at —40 mV) of WT or mutant channels after a depolarization to 4-80 mV, plotted against time. Current 
values are normalized to the current level measured before magnesium application (time 0). C, rundown time constant (i) of WT and mutant channel 
currents (n = 8-12). *p<0.05, **p<0.01, ***p<0.001 versus WT. 
doi:10.1371/journal.pone.0093255.g002 



R243H mutant as compared to WT (Figure 3B), confirming a 
decrea.sed PIP2 sensitivity. Tlie relative current amplitude after 
64 s of Ci-VSP activity was 0.35±0.04 (n= II) and 0.21 ±0.02 
(n = 6; P<0.05) for WT and R243H, respectively. When 
coexpressed with Ci-VSP, the R555C mutant channel had a very 
low current density (at time 0, I^ii = 0.16±0.08 pA/pF, p<0.001, 
n= 18) as compared to the mutant channel activity without Ci- 
VSP coexpression (I,^ii = 4.27 ± 1.64 pA/pF, n = 4). Such a Ci- 
VSP coexpression had no effect on the WT current density 
(I,^ii= 13.4±7.64 pA/pF, n = 7 and I„ii = 25.4±6.15 pA/pF, 
n=:ll without and with Ci-VSP respectively, N.S.). Consistent 
with the lower affmity of the R555C mutant to diCS-PIPa [10], 
these results suggest that Ci-VSP has enough basal activity (i.e. 
without any depolarization to -h80 mV) to decrease membrane- 
PIP2 levels and that such a basal PIPj decrease is able to affect the 
R555C channel current density without any effect on the WT 
channel. In contrast, the decrease of PIP2 levels caused by Ci-VSP 
activation had no effect on R539W channel activity (normalized I 
at 64 s = 0.99±0.06, n = 9, p<0.001 versus WT; Figure 3C). 

These three different approaches show that, contrary to R243H 
and R555C mutants, R539W is much less sensitive to variation in 
PIP2 levels than WT. A priori, the results are inconsistent with the 
similar sensitivity of R539W and R555C to diC8-PIP2 observed 
previously [10]. The impaired Mg^"*", wortmannin and Ci-VSP 
effects on R539W channel rundown questioned the validity of 
these approaches in the evaluation of PIP2 sensitivity. We wiU 



address, later in this article, why the R539W channel is less 
sensitive to PIP2 variation than WT. 

R539W is insensitive to osmolarity 

We then studied the R539W channel response to a physiological 
range of PIP2 variation, by varying extracellular osmolarity. We 
previously showed that switching from hyperosmolar to hypoos- 
molar extracellular solution leads to an increase in available 
membrane PIP2, provoking an increase in the WT KCNEl— 
KCNQ,1 current density, a shift of the activation curve towards 
negative potentials and slower deactivation kinetics (Figure 4A, 4C 
and [14]). During the same set of experiments, we measured 
R539W channel current in isoosmolar (334 mosmol.L ), hy- 
poosmolar (234 mosmol.L 70% of control osmolarity) and 
hyperosmolar extracellular solutions (434 mosmol.L" , 130% of 
control osmolarity; cf Methods). Tail current density (measured 
after a -Hl20-mV depolarization), half-activation potential and 
time constant of deactivation at —40 mV (xa^at.,) were changed 
neither by switching from isoosmolar to hypoosmolar condition, 
nor by switching from isoosmolar to hyperosmolar condition 
(Figure 4B and 4D). This lack of R539W channel osmoregulation 
is supporting the idea that the mutant is much less sensitive to 
variations in PIP2 levels than the WT channel. 

R539W is as sensitive to PKA as WT KCNEl -KCNQ1 

For some Kir channels, it has been established that phosphor- 
ylation by PKA increases their interactions with PIP2 [9,19]. 



PLOS ONE I www.plosone.org 



5 



March 2014 | Volume 9 | Issue 3 | e93255 



R539W Activity Is Regulated by Cholesterol 



0.5 s 



R539W 



0.5 s 



+80 mV 



-, -40 mV 
L -80 mV 



R555C 



Control 

64 s after Ci-VSP activation 



B 



Ci-VSP coexpressed 




1.0 



0.5 



0.0 



O 



o o 



O WT 
• R243H 

A R539W 
■ R555C 
WTw/o Ci-VSP 



O O O o O o o 



50 
time (s) 



100 



1.2 



S 0.8 



E 0.4 
o 



0.0 



\ivithogit 
Gi-Vsi= 



*** 



'5^ csi^ / 

4^ <P <P 



Figure 3. R539W is insensitive to Ci-VSP. A, representative ruptured-patch current recordings of WT or mutant channels coexpressed with the 
voltage-dependent membrane phosphatase, Ci-VSP, at the first (t = 0 s) and 9"^ (t = 64 s) step of depolarization. These currents were measured during 
the voltage-clamp protocol shown. Start-to-start interval = 8 s. The 4-80-mV depolarization also allows Ci-VSP activation. B, relative tail-current 
amplitude (measured at —40 mV) of WT or mutant channels after a depolarization to 4-80 mV, plotted against time. Current values are normalized to 
the current amplitude measured before Ci-VSP activation (time 0). C, mean relative current amplitude of WT or mutant channels measured after a 64-s 
Ci-VSP activation (n = 6-11). *p<0.05, ***p<0.001 versus WT. ^R555C already ran down before Ci-VSP activation, due to basal Ci-VSP activity, it is thus 
assimilated to 0 (n = 18). WT condition without Ci-VSP is shown in (B) and (C). 
doi:1 0.1 371 /journal.pone.0093255.g003 



Several studies have shown that KCNQ^l activity is regulated by 
PKA-dependent phosphorylation [36-38] and two studies sug- 
gested that the PKA-dependent phosphorylation of KCNQ,1 
increases its interaction with PIP2 [9,11]. If so, PKA should have 
less impact on R539W KCNEl-KGNQl current since R539W is 
less sensitive to PIP2 variation. To test this idea, cells expressing 
KCNQ^l, KCNEl and the PKA-anchoring protein yotiao, were 
exposed to a solution containing (in |j,mol/L) 400 cpt-cAMP, 10 
forskolin and 0.2 okadaic acid. Tail current density after a 1-s 
depolarization to -1-80 mV was measured in the permeabHized- 
patch configuration. cAMP had a similar effect on normalized tail 
currents of the WT and of the R539W channel (Figure 5; WT, 
1.45±0.I1 and R539W, 1.63±0.14; n= 12-13), indicating that 
the R539W mutant channel is as sensitive to PKA-dependent 
phosphorylation as the WT channel. Therefore, PKA regulation 
of KCNEl -KCNQl channel does not seem to act through a 
modulation of channel-PIPa interaction, which is consistent with 
the recent publication of Li et al, where phospho-mimetic 
mutations do not affect PIP2 dependent rundown [12]. 

Hypothesis on the molecular basis of the low sensitivity 
to PIP2 variation of R539W 

R539W and R555C mutants have fully-activated current 
amplitude similar to the WT channel [10], supporting the 
conclusion that they are similarly expressed, correctiy folded, 
and fully processed to the cell membrane. R539W and R555C 



mutants have similar sensitivity to cliC8-PIP2 on one hand [10], 
but different responses to a decrease in membrane PIP2 on the 
other hand (Figures 1, 2 and 3). Since a decrease in membrane 
PIP2 closes the WT channel [13], our data indicate that the 
R539W channel is steadily open and poorly sensitive to membrane 
PIP2 variation. Thus, R539W open state is not or less dependent 
on a direct or indirect interaction with PIP2, as compared to WT. 

To better understand the mechanism by which R539W is 
desensitized to variations in PIP2 levels, we used molecular 
modeling to analyze the structure of the WT and mutant 
Q.QARKPYDV-R/W539-DVIEQYSQG fragments of the 
KCNQ^l C-terminus and their potential interaction with the 
membrane bUayer. In the absence of a fuU 3D model of the 
channel, only partial models were used to reach a working 
hypothesis. 

Three dimensional structures of the fragments were calculated 
with PepLook [30]. For both fragments the 99 lower energy 
models were sorted and compared. For the WT sequence, the 99 
models had a similar conformation i.e. a fi-extended polar hairpin. 
R539 was located in the U-turn (Figure 6B). Models were clustered 
into three leads on the basis of RMSd less than 1 A (Figure 6A), 
Leads were stabilized by aromatic (Y536/Y545) side chain 
interactions or by NH...Tt and NH...O H-bonds but mainly by 
side chain (R533-E543 and/or K534-D540) salt bridges. Of note, 
mutations of R533 and E543 are associated with LQTl [39], 
consistent with their role in the hairpin structure stabilization. The 
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Figure 4. R539W is insensitive to osmoiarity. A, B, superimposed representative permeabilized-patch recordings of WT (A) and R539W (B) 
KCNE1-KCNQ1 concatemer currents, respectively, measured in hyper-, iso- and hypoosmotic conditions using the voltage protocol shown in the 
insert. C, D, averaged tail-current density (in pA/pF), V0.5 (in mV) and Tdeact (in ms) measured at —40 mV after a depolarization step to 4-80 mV for WT 
(C) and 4-1 20 mV for R539W (D) channel, in hyperosmotic (hyper), control (iso), and hypoosmotic solutions (hypo) (n = 1 0); same voltage protocol as in 
A. *p<0.05, ***p<0.001 (one-way ANOVA for repeated measures). Protocol and experimental conditions are similar to those used in our previous 
study analyzing the osmoregulation of KCNE1-KCNQ1 [14]. 
doi:1 0.1 371/journal.pone.0093255.g004 




°0.2s 



Figure 5. R539W is sensitive to PKA. A, representative permeabilized-patch current recordings of WT or R539W channels (co-transfected with 
yotiao) measured during the voltage-clamp protocol shown. B, mean time course of channel activation by 400 cpt-cAiVlP, 10 |iM forskolin and 
0.2 |iM okadaic acid (cAMP) for WT or R539W KCNE1-KCNQ1 concatemer channel tail currents measured at -40 mV after a 1 s depolarization to 
4-80 mV and normalized to the current value before cAMP application (n = 12-13). 
doi:1 0.1 371/journal.pone.0093255.g005 
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lead models were polar structures with a mean hydrophobic to 
hydrophilic accessible surface ratio of 0.73±0.20 (Figure 6C). 
Because of its hydrophilicity, the fragment should map to the 
channel protein surface rather than in its core, being accessible to 
external partnership. In addition, due to its high content in 
positively-charged residues (R533, K534 and R539) and especially 
to the fact that R539 was accessible in all models, the fragment 
might be able to interact with negative charges of membrane 
phospholipids. As a basis for this hypothesis, we found that the best 
position of the WT fragment models was to be adsorbed on the 
membrane surface, with R539 diving in the interface (Figure 7A). 
In that position, what could be the R539 membrane partner? In 
Figure 7D, we show the optimal positions of PIP2, cholesterol and 
DOPC in the same membrane slab. The depth of insertion of the 
WT fragment models argues for a possible interaction of R539 
with the C = O moiety of PIP2 (Figure 7D-7F). 

For the R539W mutant, changing R to W changed the 3D 
structure but not the extended conformation (Figure 6A compared 
with 6D). Four lead models were sorted (Figure 6D) and W539, 
like R539, was a structural protuberance. The hydrophobicity 
profiles of the WT and R539W models demonstrate that changing 
R to W decreased the hydrophilicity (Figure 6F to be compared 
with Figures 6C, 7B and 7C). The R539W mutant had a mean 
hydrophobic to hydrophilic accessible surface ratio of 0.81+0.10 
supporting the conclusion that the fragment was still polar and that 
the R539 polar protuberance of the native fragment was now a 
W539 hydrophobic protuberance. The best position of the 
R539W mutant channel fragment in the membrane /water slab 
was the membrane interface similar to the WT fragment 
(Figures 7A-7C). Even though fragment positions are similar, 
their partners in the membrane should be different. The WT 
fragment could have a polar R-phospholipid interaction whereas 
the mutant fragment would more likely favor an interaction of the 
apolar W539 with a more hydrophobic partner. 

Among natural membrane lipids, cholesterol has been impli- 
cated in direct interactions with some channels and in the 
regulation of their activities [40,41]. Various studies have 



suggested that protein-cholesterol interactions implicate trypto- 
phan [42,43]. Interestingly, in another ion channel, a mutation of 
a cysteine to a tryptophan in a lipid-exposed position enhances the 
effect of cholesterol [42] . Tryptophan can interact with cholesterol 
in a rigid cycle-to-cycle hydrophobic fit stabilized by OH 
(cholesterol) - 7i (tryptophan cycle) or by NH (tryptophan cycle) - 
O (cholesterol) interactions if the two partners are at the same 
depth of insertion in the membrane. Our calculations indicate that 
the cholesterol polar head and tryptophan are at the same level in 
the membrane and thus justify the hypothesis that changing R539 
to W could switch a polar R-PIP2 interaction to a W-cholesterol 
interaction. 

R539W interacts with membrane cholesterol 

To test this hypothesis, we probed whether the decreased 
sensitivity of the R539W mutant to PIP2 was associated with a gain 
in cholesterol sensitivity. We tried three diflferent approaches: first, 
we measured the Mg '^"-induced current rundown of the R539F 
KCNE 1 -KCNQ_1 channel, phenylalanine being less hydrophobic 
and bulky than tryptophan and thus likely to be a less efficient 
membrane anchor than tryptophan. Second, we compared the 
effect of depleting the membrane cholesterol by 2-hydroxypropyl- 
fi-cyclodextrin (cyclodextrin) on WT and R539W channel current 
[44] . Since cyclodextrin is not very specific, we used another way 
of depleting membrane cholesterol. We selected triparanol 
because it has been shown to decrease cholesterol in a similar 
model [6]. In all these approaches we measured the relative tail- 
current amplitude of WT and mutant channels after a depolar- 
ization to +80 mV, during 1.1-mmol/L free Mg^"*" application in 
the giant-patch configuration. 

As expected, the Mg "^-induced current rundown of the R539F 
mutant channel was faster than the rundown of the R539W 
mutant channel activity and as fast as the rundown of WT current 
(Figures 8A and 8B). This behavior is consistent with a lower 
anchoring of the F539 on the membrane as compared to that of 
W539. The fact that R539F mutant channel is not running down 



W539 



Figure 6. PepLook models of the 19-aa sequence surrounding R539 and W539. The sequences QQARKPYDVR539DVIEQYSQG and 
QQARKPYDVW539DVIEQYSQG were used to calculate amphipathic 3D structure models in water. A, D, the 99 PepLook models of low energy were 
clustered into three WT (A) and four R539W (D) lead models on the basis of backbone RIVlSd<1 A. Ribbon structures of these lead models are fitted in 
(A and D) demonstrating their large similarity. B, E best models (Prime), with R539 in CPK for the WT fragment (B) and W539 in CPK for the mutant 
fragment (E). Both residues are protruding in a pin loop-like structure. C, F, hydrophobicity profiles of the WT and R539W Primes are visualized by the 
hydrophobic (brown) and hydrophilic (green) isopotential surface (+0.1 kcal/mol) around the molecule. The MHP profiles demonstrate, first, the 
rather important hydrophilicity of the models, and second, that changing R to W, transforms a very polar protuberance to an apolar one. 
doi:1 0.1 371/journal.pone.0093255.g006 
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Figure 7. Relative position of a series of molecules in the 
membrane. All molecules are shown at their optimal position in the 
IMPALA slab after they were systematically tested at every A across a 
water/membrane continuous layer. In all plots, the yellow grid 
represents the membrane center; the purple grid, the averaged lipid 
head/acyl chain interface, and the pink grid the averaged lipid/water 
interface. A, PepLook model of the WT (left) and R539W (right) peptides. 
In both cases, residue 539 is imbedded in the membrane. B, MHP 
(Molecular Hydrophobicity Potential) profiles of fragments. R539 is 
responsible for a large hydrophilic (green) protuberance, whereas W539 
is a hydrophobic (brown) protuberance in the membrane. C, IVIEP 
(molecular electrostatic profile) of the same molecules showing the e- 
attractivity of the R539 protuberance. D, E and F, same as in A, B and C 
except that molecules are, from left to right, an extended fluid form of 
PIP2, cholesterol, and a fluid form of DOPC. They highlight the relative 
position of PIP2 and peptides and the fact that the more hydrophobic 
cholesterol is embedded in the membranes. 
doi:1 0.1 371 /journal.pone.0093255.g007 

faster than WT suggests that phenylalanine still interacts with 
cholesterol although less than the tryptophan does. 

Second, we compared the Mg^"^-induced rundown of the 
R539W mutant current without and with 2 mmol/L of cyclodex- 
trin. As shown in Figures 8C and 8D, cyclodextrin pre-treatment 
accelerated the R539W rundown. Thus, when membrane 
cholesterol is decreased, the R539W channel is more sensitive to 
PIP2 variation. In WT channels, only activation kinetics were 
affected by cholesterol depletion, as shown in Figure S 1 in File SI 
and consistent with a previous work [6], but rundown was not 
affected (Figure 8C and 8D). As opposed to WT channels, R539W 
activation kinetics were not modified by cyclodextrin (Figure S 1 in 
File SI). 

Finally, we compared the Mg^'''-induced rundown of the 
R539W mutant current without and with 10 |J,mol/L of tripar- 
anol. Like cyclodextrin, triparanol pre-treatment accelerated the 
R539W rundown without any significant change in the WT 
rundown (Figure 8E and 8F). These results confirm that the 
R539W channel is more sensitive to PIP2 variation when 
membrane cholesterol is decreased. In addition, as for cyclodex- 
trin, in WT channels, only activation kinetics were affected by the 
triparanol pre-treatment. In R539W mutant channels, activation 
kinetics were unchanged by the triparanol pre-treatment (Figure 



S2 in File SI), consistent with the results using cyclodextrin to 
deplete membrane cholesterol. These results reinforce the idea 
that interaction of cholesterol with WT and R539W channels is 
quite different. 

Altogether, data from experiments manipulating the residue at 
position 539, and cholesterol levels in the membrane, are in 
agreement with the hypothesis that the tryptophan in position 539 
interacts (or induces an interaction) with membrane cholesterol. 
This channel-cholesterol interaction might overrule the channel- 
PIP2 interaction to stabilize the channel open-state. 

Discussion 

In previous studies, we have shown that KCNEl-KCNQ^l 
channel open-state is stabilized by PIP2 and impairment of this 
stabilization by arginine neutralization at position 243, 539 or 555 
in KCNQ^l is correlated with the long QT syndrome [10,13]. For 
the three mutants, higher diC8-PIP2 concentrations than for WT 
were needed to stabilize the open state after the channel activity 
had run down, suggesting a decrease in interaction with PIP2 [10]. 
The R539W and R555C mutations are localized in the cytosolic 
C-terminus (CTD) [45]. The supposed interaction of arginines 539 
and 555 with PIP2 suggests that they are situated on the 
membrane-cytosol interface, which may be surprising since they 
are located in the middle of the distal half of the KCNQl cytosolic 
CTD [45]. In a recent crystallographic study, Hansen et al. 
described that PIP2 mediates docking of the whole CTD to the 
transmembrane module and subsequent opening of the inner helix 
gate of the Kir2.2 channel [46] . Thereby, the KCNQ,1 distal CTD 
might come close in order to interact with the membrane, via 
interactions such as R539-PIP2 and R555-PIP2, allowing the CTD 
to be in the vicinity of the pore domain for modulating its opening. 
However, further crystallographic studies should address this 
hypothesis. 

In the present study, the most striking result was that decreasing 
endogenous PIP2 has a very small effect on the open-state stability 
of the R539W mutant. R539W is paradoxically less sensitive to a 
decrease in PIP2 than WT. This is in contrast with many channel 
mutations that reduce their afiinity for PIP2 or diC8-PIP2 [47,48], 
including the two other KCNQl mutations R243H and R555C. 
Here we suggest that the paradoxical behavior of R539W is due to 
the stabilizing effect of tryptophan-cholesterol interaction that 
replaces the arginine-PIP2 interaction. Indeed, severe cholesterol 
depletion of the membrane restored the PIP2 sensitivity of the 
channel. After cyclodextrin or triparanol pre-treatment, R539W 
rundown kinetics during Mg^"*" application were restored to values 
similar to WT, suggesting that membrane cholesterol depletion 
abolished the R539W open pore stabilization, now only 
maintaining its open stabilization through other PIP2-interactiiig 
residues (such as R243 and R555). It is important to note that - 
when expressed in Xenopus oocytes - R539W behaved exactly as 
R555C [12]. In that model, both excision-induced rundown and 
PIP2 sensitivity were the same for R555C and R539W. This 
apparent inconsistency between models may be due to different 
cholesterol levels around the channel in X. oocytes versus COS-7 
cells, with higher levels in the latter, slowing down R539W 
rundown. 

For Kir4.1, Hibino and Kurachi showed that cyclodextrin 
abolishes the function of Kir4.1 channels in HEK293 cells [49]. 
One of the proposed mechanistic hypotheses suggests that 
cholesterol influences the conformation of Kir4. 1 and activates it 
by specific protein-cholesterol interactions [49] . This hypothesis is 
supported by some studies performed on other channels. Singh et 
al. demonstrated a direct effect of cholesterol on the activity of 
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Figure 8. Effect of intracellular magnesium on WT, R539W and R539F mutant channels, and after membrane cholesterol depletion. 

A, C, E, relative tail-current amplitude (at —40 mV) of WT or mutant channels measured after a depolarization to +80 mV plotted against time during 
a 1 .1 mmol/L free Mg^* application on a giant-patch. Start-to-start interval = 2 s. Current values are normalized to the current level measured before 
magnesium application (time 0). Cholesterol depletion was induced by 1 hour of 2 mmol/L cyclodextrin (cycio) or 24 hours of 10 |imol/L triparanol 
(tripa) pre-treatment. B, D, F, mean rundown time constant (i) of WT and mutant channels (B), and with and without 2 mmol/L cyclodextrin (D, n = 9- 
13) or 10 |imol/L triparanol (F, n = 9-15). *p<0.05, **p<0.01 versus WT. DMSO in which triparanol was diluted (E and F) has an effect on rundown 
kinetics (t = 15.4±1.6, n = 7 versus t = 35.0±3.65, n = 15, without and with DIVISO pre-treatment respectively, p<0.01). 
doi:10.1371/journal.pone.0093255.g008 



KirBac 1 . 1 [40] . They showed that changes in membrane fluidity 
could not account for the effects of the sterols on KirBac 1.1 
activity, and their data strongly pointed to a direct cholesterol- 
channel interaction. The present study is also supporting a direct 
and specific cholesterol-channel interaction. 

We have previously reported that KCNE1-KCNQ_1 activity is 
regulated by intracellular ATP in addition to PIPj [13]. More 
precisely, Li et al. recently showed that ATP binds to the cytosolic 
domain and promotes pore opening [50]. Since the patch-clamp 
experiments were performed without ATP in the intracellular 
solution, the gradual decrease of channel activity is probably the 



sum of the PIP2- and the ATP-dependent channel rundown. This 
is especially true in giant-patch experiments, during which ATP is 
diffusing out of the patch very quickly. Even though the regulation 
of KCNQl by PIP2 and by ATP concerns two distinct 
mechanisms [13], and distinct binding sites [51], there may be 
some interplay between both types of regulation. For instance, it 
has been shown that variations in intracellular ATP provoke 
variations in membrane PIP2 levels [25,52]. In that context, using 
R539W - which is much less PIP2-sensitive - may reveal to what 
extent PIP2-dependent rundown influences the ATP-dependent 
rundown. 
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KCNC)1 and KCNEl localize in lipid rafts [53-55] as well as 
the G-protein-dependent machinery involved in |3-adrenergic 
signaUng [56] . From the literature, it appears that cholesterol may 
be involved in KCNQ,1 expression levels and channel-complex 
trafficking [5,57] and gating [6]. Here we show that cholesterol 
does not stabilize WT KCNE 1 -KCNQl channel opening per se. 
Substitution of R539 by W introduces a new sensitivity to acute 
cholesterol changes, possibly by favoring a direct binding of this 
residue to cholesterol, which is abundant in lipid rafts. In Kir2.1 
channels, Epshtein et al. identified a specific region that plays a 
critical role in the sensitivity of these channels to cholesterol [58]. 
This region is in the CTD, consistent with the location of R539 in 
KCNCil. Nonetheless, PIPj still regulates partially the R539W 
mutant, since in the presence of cholesterol, severe depletion of 
PIP2 eventually induced open-state destabUization. This indicates 
that other PIP2 interacting residues (such as R243 and R555) are 
still needed for channel opening. This severe depletion of PIP2 
allowed unmasking of the reduced diC8-PIP2 sensiti\ity of the 
R539W mutant, since higher doses of diC8-PIP2 are necessary to 
activate the channel [10]. 

We inferred the PIPj sensitivity of WT and mutated channels 
from the measurement of the rundown kinetics of the correspond- 
ing current during a decrease in available PIP2.The decrease in 
PIP2 was triggered using three common approaches: extracellular 
wortmannin application on permeabilized-patch configuration, 
intracellular Mg^"*" application on excised-patch configuration, and 
Ci-VSP coexpression on ruptured-patch configuration 
[10,18,27,28,34]. Although the R539W mutant showed a 
decreased affinity to diC8-PIP2 compared to WT [10], all these 
three approaches show that this mutant is less sensitive to a PIP2 
decrease, which would classically suggest a higher aflOnity to PIP2. 
This peculiar behavior suggests that caution should be taken when 
using the current rundown in the evaluation of a channel PIP2 
affinity, because, in some cases (e.g. R539W), the mutation that 
disrupts the interaction with PIP2 may stabilize the channel in the 
open state through another mechanism. 

In conclusion, the present study shows that a channel mutation 
can induce or favor an interaction between the channel and a 
membrane component. We speculate that this PIP2-independent 
tryptophan-cholesterol interaction is responsible for the difference 
in phenotypes observed for patients with long QT syndrome 
caused by R539W and R555C mutations of KCNQ.1. Until now. 
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